Dissipative processes in superfluid neutron stars 
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Abstract. We present some results about a novel damping mechanism of r-mode oscillations in neutron stars due to processes 
that change the number of protons, neutrons and electrons. Deviations from equilibrium of the number densities of the various 
species lead to the appearance in the Euler equations of the system of a dissipative mechanism, the so-called rocket effect. The 
evolution of the r-mode oscillations of a rotating neutron star are influenced by the rocket effect and we present estimates of 
the corresponding damping timescales. In the description of the system we employ a two-fluid model, with one fluid consisting 
of all the charged components locked together by the electromagnetic interaction, while the second fluid consists of superfluid 
neutrons. Both components can oscillate however the rocket effect can only efficiently damp the countermoving r-mode 
oscillations, with the two fluids oscillating out of phase. In our analysis we include the mutual friction dissipative process 
between the neutron superfluid and the charged component. We neglect the interaction between the two r-mode oscillations 
as well as effects related with the crust of the star. Moreover, we use a simplified model of neutron star assuming a uniform 
mass distribution. 
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INTRODUCTION 

Fermionic superfluidity might be realized in the interior 
of neutron stars where the temperature is low and the 
typical energy scale of particles is extremely high. In 
particular, in the inner crust of standard neutron stars the 
attractive interaction between neutrons can lead to the 
formation of a BCS condensate. 

A way to detect or to discard the presence of a super- 
fluid phase in compact stars consists in studying the evo- 
lution of the r-mode oscillations [1]. R-modes are non- 
radial oscillations of a fluid with the Coriolis force act- 
ing as the restoring force. They provide a severe limita- 
tion on the rotation frequency of a star through coupling 
to gravitational radiation. When dissipative phenomena 
damp these oscillations the star can rotate without losing 
angular momentum to gravitational radiation. If dissipa- 
tive phenomena are not strong enough, the r-mode oscil- 
lations will grow exponentially fast in time and the star 
will keep slowing down until some dissipation mecha- 
nism is able to damp the r-mode oscillations. In this way 
one can put some constrains on the stellar structure rul- 
ing out phases that do not have large enough viscosity. 
For such studies it is necessary to consider in detail all 
the dissipative processes and the interactions among the 
various layers of a star 

In Ref. [2] we have studied a novel dissipative pro- 
cess associated with the change in the number of pro- 
tons, neutrons and electrons. In real neutron stars these 



processes can take place in the outer core and in the in- 
ner crust of the star and are related to beta decays and in- 
teractions between the neutron fluid and the crust. Both 
these processes lead to the appearance of a dissipative 
force (the so-called rocket term) in the Euler equations 
of the system. This force is due to the fact that when two 
or more fluids move with different velocities a change 
of one component into the other results in a momentum 
transfer between the fluids. This change in momentum 
is not reversible, because it is always the faster moving 
fluid that will lose momentum. The resulting dissipative 
force is proportional to the mass rate change, and to the 
relative velocity between the fluids. The name "rocket ef- 
fect" reminds that the same phenomenon takes place in 
the dynamical evolution of a rocket whose mass is chang- 
ing in time as it consumes its fuel. 

As far as we know, the dissipative force due to the 
rocket term has not been considered in the context of r- 
mode oscillations. In the hydrodynamical equations cor- 
responding to the mass conservation laws, it is in general 
assumed that the neutron, proton and electron compo- 
nents are separately conserved quantities. Assuming that 
the change in the particle densities is due to out of equi- 
librium direct Urea processes, we have determined the 
typical timescale associated with the rocket effect and 
we have found that it is sufficiently short to damp coun- 
termoving r-mode oscillations, with the normal and su- 
perfluid component oscillating out of phase. 



ROCKET EFFECT 

For a system consisting of neutrons, protons and elec- 
trons, the mass conservation law is given by, see e.g. [3], 

dtpx + yi{pxv'J = r,, (1) 

where is the particle mass creation rate per unit vol- 
ume and the index x = n,p,e refer to the particle species, 
that is, neutrons, protons and electrons. In these equa- 
tions we have considered that some process can convert 
neutrons in protons and electrons and vice versa. There- 
fore, we are assuming that the various components are 
not separately conserved. 

One possible mechanism leading to a change in the 
particle number densities is the crust-core transfusion, a 
process that may take place when the ionic constituent of 
the crust are squeezed by the underlying superfluid and 
part of their hadronic content is released and augments 
the superfluid component. The opposite mechanism, re- 
lated to a reduction of the pressure leads to the capture of 
protons and neutrons by the ions of the crust. 

A different mechanism leading to a change in the 
particle number densities is the direct Urea process 

n^p + e~ + Ve, p + e~ ^n + Ve- (2) 

It was found in Ref. [4], that for certain realistic equa- 
tions of state the direct Urea processes are allowed when 
the star density exceeds the nuclear saturation density 
po = 2.8 X 10^^ g cm^^, and the proton fraction exceeds 
the threshold value xfp = ^. We shall restrict to consider 
such processes which might be realized in the interior of 
massive neutron stars. 

The three particle creation rates are not independent 
quantities, because charge and baryon number conserva- 
tions imply that = = — F„, meaning that only one 
creation rate is independent. 

It is possible to simplify the treatment of the sys- 
tem considering that electrons and protons are locked 
together by the electromagnetic interaction. Thus they 
move with the same velocity. Moreover, charge neutral- 
ity implies that ne=np. Therefore, electrons and protons 
can be treated as a single charge-neutral fluid and hence- 
forth we shall refer to this fluid as the "charged" compo- 
nent. The two fluids have mass densities 

Pn = nin rin and Pc = m„ tic , (3) 

where m„ = mp+ nig and n^ = np = ng. The Euler equa- 
tions obeyed by the two fluids are given by 

(c'* + v^V,.)K + e„w,-) + V,.(/l„ + 4>) + e„w^V,.v^ = 0, 
{dt+vi:Vj){v'i-ecWi) + Vi{fic + ^)-ecwjViv'j = Rt 

where <I> is the gravitational potential, i, j label the space 
components; we have defined a chemical potential by 



mass p.x = IJ-x/m„, and w = Vc — v„ represents the relative 
velocity between the two fluids. The quantities e„ and Ec 
are the entrainment parameters, that are related to the fact 
that momenta and velocities of quasiparticles may not be 
aligned [5]. 

The term on the right hand side of Eq. (4) is given by 
RT = {l-en-ec)—Wi (4) 

Pc 

and represents the rocket term. In the analysis of the pos- 
sible dissipative mechanisms of star oscillations this term 
is usually neglected. Indeed, it is in general assumed that 
the neutron, proton and electron numbers are separately 
conserved quantities, that is Fp = Fg = F„ = 0. In Ref. [2] 
we have considered the effect of this term and obtained 
the corresponding damping timescale, assuming that the 
change in the number densities is due to direct Urea pro- 
cesses. 



Stability analysis 

In our analysis we have considered two different r- 
mode oscillations. One is associated predominantly with 
toroidal comoving displacements of the normal and su- 
perfluid components of the system. The second type of 
r-mode oscillation is dominated by toroidal countermov- 
ing (out of phase) displacements of the two components. 
We refer to these oscillations as "standard" r-modes and 
as "superfluid" r-modes, respectively. These two modes 
decouple for a star made by uniform and incompressible 
matter, and we shall restrict to treat such a case. 

As in Ref. [6], we have studied the linearized hydrody- 
namical equations for the perturbations around an equi- 
librium configuration of a neutron star rotating with con- 
stant angular velocity il, and we assume that the back- 
ground configuration is such that the two fluids move 
with the same velocity, thus at equilibrium w = 0. 

For simpUcity we have considered a toy-model neu- 
tron star with uniform density p = 2.5po, with radius 
R=10 km and mass M ~ 1 .47 Mq, and no crust. 

In Fig. 1 we report the damping timescale associated 
to the rocket effect as a function of the entrainment 
parameter e, see Ref. [2], for standard r-modes. The 
typical timescale of gravitational wave emission is of the 
order of tens of seconds, therefore the rocket effect is not 
able to efficiently damp these oscillations. 

In Fig. 2 we report the results for the superfluid r- 
mode "instability window" for a star with uniform den- 
sity, p = 2.5po, and R = \0 km. It is obtained compar- 
ing the various dissipative timescales with the growth 
timescale of gravitational wave radiation. Stars above the 
various lines are unstable. The dashed red Une represents 
the instability window in the absence of the rocket effect. 
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FIGURE 1. Damping time associated with the rocket term 
decreases with increasing entrainment parameter e, see Ref . [2] 
for more details. We have taken T = lO^K and a critical super- 
fluid temperature Tc = lO'^K. 



see e.g. [7], for two different values of the entrainment 
parameter. 
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FIGURE 2. Instability window of the superfluid r-modes of 
a star with uniform density p = 2.5po, with radius R= \Q km 
and mass M ~ 1.47M0. The dashed red line represents the 
instability window in the absence of the rocket term. The full 
blue line represents the instability window with the inclusion 
of the rocket term. The horizontal full lines correspond to the 
effect of the mutual friction for £ = 0.0002, lower line, and 
e = 0.002 upper line. In our simplified model of star, the mutual 
friction is independent of the temperature. 
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The region above the dashed red line is unstable when 
only shear and bulk viscosity damping are considered. 
At low temperature, shear viscosity is the dominant dis- 
sipative mechanism. With increasing temperature shear 
viscosity is less efficient and it can damp r-mode oscilla- 
tions for smaller and smaller values of the frequency. The 
behavior of the bulk viscosity is the opposite and starts 
to damp r-mode oscillations for temperatures of the order 
of lO'" K. The full blue line represents the effect of the 
rocket term. It has a behavior qualitatively similar to the 
bulk viscosity, but it becomes effective at smaller tem- 
peratures. Therefore, the instability window with the in- 
clusion of the rocket term is much reduced. The full hor- 
izontal lines correspond to the effect of mutual friction. 



